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Abstract 

By generalizing the auxiliary field term in the Lagrangian of simplicial chiral models 
on a (d — 1) -dimensional simplex, the generalized simplicial chiral models has been 
introduced in jl]]. These models can be solved analytically only in d = and d = 2 
cases at large-N limit. In d = case, we calculate the eigenvalue density function 
in strong regime and show that the partition function computed from this density 
function is consistent with one calculated by path integration directly. In d = 2 case, 
it is shown that all V = Tr(AA T ) n models have a third order phase transition, same as 
the 2-dimensional Yang-Mills theory. 
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1 Introduction 



The 1/N expansion of matrix-valued field theories is probably the most important nonper- 
turbative and nonnumerical theoretical tool presently available in the study of such models 
as non-Abelian gauge theories and two-dimensional quantum gravity. The first major result 
concerning the large-N limit of matrix theories was due to 't Hooft, who made the crucial 
observation that in the 1/N expansion of continuum gauge theories, the set of Feynman 
diagrams contributing to any given order admits a simple topological interpretation ||. 

Unfortunately, our knowledge of exact solution for the large-N limit is limited to a small 
number of few-matrix systems, and it becomes even smaller if we restrict ourselves to the 
case of unitary matrix fields, which is relevant to the problem of lattice QCD. After the 
solution of Gross and Witten || of the single link problem and its generalizations, exact 
results were obtained only for the external field problem [||, and a few toy models (L=3,4 
chiral chains) [^, |7[]. Therefore extending the number of solved few-matrix systems is an 
important progress from different points of view: first from pure theoretical reason. Second 
by noticing that any few-matrix system have a reinterpretation, via double scaling limit, 
as some different kind of matter coupled to two-dimensional quantum gravity. And third 
because every few-matrix system involving unitary matrices can be reinterpreted as the 
generating functional for a class of integrals over unitary groups, and these integrals in turn 
are the essential missing ingredient in the context of a complete algorithmization of the 
strong coupling expansion of many interesting models ||. 

One of the interesting classes of the few-matrix models is the simplicial chiral model 
(SCM) 0. In this model, each unitary matrix interacts in a fully symmetric way with all other 
matrices, by preserving the global chiral invariance. The resulting system can be described as 
a chiral model on (d — l)-dimensional simplex. A simplex is formed by connecting d vertices 
by d(d — 1)/2 links. In ref.[10|, where the large-N saddle-point equations for density function 
p(z) have been found, the authors have introduced a single auxiliary variable A (a complex 
matrix) to decouple the matrix interaction, and by performing the single-link external field 
integral, the saddle-point equation in strong and weak regions have been found. In d — 2, 
d = 4, and d — > oo, the saddle-point equation has been solved analytically and the phase 
transition of the model has been specified. In the action of SCM in terms of the matrix A, 
this matrix appears as a Tr(AA^) term. 

On the other hand, it is well known that the pure 2-dimensional Yang-Mills theory 
(YM2) can be represented by the Lagrangian iTi(BF) + Ty(B 2 ), in which F is the usual 
field-strength tensor and B is an auxiliary pseudo-scalar field. Many of the main properties 
of YM 2 does not change if one considers the generalized two-dimensional YM 2 (gYM 2 ), 
which can be found by replacing Tr(i? 2 ) term in YM 2 action by an arbitrary class function 
f(B) |~T|. The partition function of these theories have been fully studied in different contexts 



1 



in refs.|jLT|-[|i~3| and the phase structure of some of the specific examples has been studied in 
refs. |14| and [15| . In all the studied cases, it is seen that the models have third-order phase 
transition, which is the same behavior as YM 2 . The crucial point in this area is that such 
investigation for continuum gYM2 is very complicated and there is no general result for the 
phase transition of an arbitrary gYM 2 theory. 

In ref. the procedure used to obtain the gYM2 from YM2, i.e. Tr(I? 2 ) — > f(B), has 
been used in SCM to introduce the generalized simplicial chiral model (gSCM). That is, 
the Tr (AA') term in SCM has been replaced by an arbitrary class function V(AA'). The 
large-N limit of the model has been discussed and the saddle-point equations in strong and 
weak regimes have been found. Note that as the SCM at d = 2 is the discrete version of 
YM 2 [p!0| |) we can consider the gSCM, at d = 2, as the equivalent matrix theory of gYM 2 . 
The phase structure of the model for d = 2 and V = Tr(AA^) n (n = 2, 3, 4) is also obtained 
in|fl. 

In this paper we want to complete our investigation of d = 2 gSCM and show that 
all V = Tr(AA^) n models (with n an arbitrary positive integer) have third-order phase 
transition. We think that it is an important result as it indirectly proves the equivalence of 
all Ti(B n ) gYM 2 theories with YM 2 , from the phase structure point of view. The plane of 
the paper is as follows. In section 2, we bring a brief review of SCM and gSCM and also 
the saddle-point equations in weak and strong regimes. In section 3, we focus on d = 
case and show that the calculation of the partition function of the model, described in 
terms of A fields, leads to the trivial result which obtained from the main action. This is a 
consistency check of our procedure. Finally in section 4 we investigate the phase structure 
of V = Tt(AA^) u gSCM in d = 2 and show that all of these models have third-order phase 
transition. We also discuss the variation of the discontinuity of f3 2 dC(f3, N)/df3 with n at 
the critical point (C is the specific heat and (3 = (g 2 N)~ l , where g is the coupling constant). 



2 The gSCM 

If we assign a U(iV) matrix to each vertex of a (d—1) -dimensional simplex, then the partition 
function of simplicial chiral models is defined by || 

Z d (/3,N) = J fldU.expiNPY, E MUiUj + UjUj)}. (1) 

i=l i=l j=i+l 

To find this partition function, the authors of || have introduced a single auxiliary variable 
A to decouple the matrix interactions. The resulting partition function is 

Zd = Z d / Z , (2) 
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where 

/d 
J] dUidAexp{-NpTrAA^ + N/3TrA ^ U} + NfiTrA^ J2 U i~ N 2 fid}. (3) 

i=l i i 

Performing the single-link integral over the matrices Ui 

e NW(BBt) = J dUexp [ NTr ^U + IJ\ B ^ ( 4 ) 

we obtain 

Z d = J dAexp{-N(3TrAA^ + NdW((3 2 AA^) - N 2 (3d}. (5) 
The gSCM is defined through the action U 

Zd,v = Zdy/Z y , (6) 

in which 

Z dy = J dAexp{-N(3V{AA^) + NdW{(3 2 AA^) - N 2 (3d}, (7) 
where ^(AA^) is an arbitrary class function of AA': 

V(AA^) = a n Tr(AA*) n . (8) 

n=l 

In special case a n = S nt x, the gSCM reduces to SCM. If we denote the eigenvalues of the 
Hermitian semi-positive-definite matrix /3 2 AA^ by Xj's, and at the large-N limit, which we 
are interested in, using the expression W(xi) in the strong and weak regimes |3[ it can be 
shown that the saddle-point equation for the eigenvalue density function p(z) (zi = + c) 
is!: 

In the above equation V indicates the principal value of integral, and the parameters a and 
b must be determined dynamically. The normalization condition of p(z) is 

b 

p{z')dz = 1. (10) 
In the weak coupling regime (/3 > /3 C ), c w = and the density function p(z) must satisfy 

2. (11) 

J a Z 

In the strong coupling regime (/3 < /5 C ), c s = a 2 and p(z) must satisfy 

V# = 2. (12) 

a Z 
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Using the standard method of solving the integral equation [16|, it can be shown that the 
density function p(z) in weak regime satisfies the following equation 



sj{b-z){z-a) 



^ I m,p,q P 

d-2 f b dy p w (y) 
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f'JV_ P.W for ,8 > A, (13) 

!/ + ^(6 + t/)(t/ + o.) J 



where C n = (2n — l)!!/(2 n n!), and in strong regime satisfies 



/m^i for/3<(3c , (14) 

2 J a y + z\ y + b y J 

where _B m = (2m — 3)!!/(2 m m!) (i?o = —1) IJ- Also by investigating the behavior of the 
integrals at z — > oo, it can be shown that the parameters a and b m. f3 > (3 C must be 
determined from the following equations 

/02„_l CrnQ'ln-m-\Q m b 2n " 1 — 2 = 0, (15) 

and 

E ^C m C 2n . m a m ^ m - (a + 6) = 1, (16) 

n,m P 

and in strong regime, (3 < /3 C , from the equations 

n,m,p y \ f 

E |^T(- a2 ) P f n " ^Can-^a^fi**- 5 ^ + a - 6 = 1. (18) 

i,m,p P \ P J 

Finally if we denote the internal energy per unit link by U d y, we have by definition 

d -^U dy = ^ 2 ^Z dy ^N). (19) 

After some calculation, it can be shown that at large-N limit, the internal energies in weak 
and strong regimes for V = Tr(AA^) n are: 

, . On — 1 r b 1 1 

d(d - 1)US = jf z 2n p w {z)dz - d +(--2)-, (20) 

and 

, . On — 1 r b 1 1 

d(d - l)U { d % = J a (z 2 - a 2 r Ps (z)dz - d + (- - 2)i (21) 

These equations can be used to deduce the order of phase transition of the models. 



E 7^T(-« 2 ) P r 1 )B m C 2n - 2p - m -ia m b 2n - 2 ^- 1 + 2 = 0, (17) 



and 
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3 d = gSCM with V = Tr(AAt) 



It is clear from the definition of Z^v m e Q-@ that Z 0i y = 1, however Z y is not trivial. To 
see this, let us focus on V = Tr(yL4^) n from now on. It is not difficult to show that Zo >n 
from eq.(^) is H 

Zo,n ~ exp [— ^ • (22) 

Now to check the procedure introduced in the last section for studying the theory in the 
large-N limit, it is instructive to reproduce this result by using the density function p(z) at 
d = 0. 

At d = 0, the saddle-point equation @ in weak regime (c = 0), for the case = 8k >n , 
becomes 

w^^ =v i¥=^ ij • " >ft (23 > 

and in the strong regime (c = a 2 ) becomes 

n <z 2 - a 2 )- 1 = P f b ^\dz' , /?<&. (24) 



2p 2n - lK ' J a Z 2 -Z' 2 

Let us focus on (3 < (3 C case. To make the integral equation fl2~4] ) more conventional, we use 
the change of variable z — > A = z 2 , with density function p s (A) satisfies 

p s (X)dX = p s (z')dz'. (25) 

So 

^zr(A-a 2 r- = p/J^l^ , (26) 
Now consider the function H S (X) in complex- A plane 

H.(\) = f ^-M. (27) 



la? \~i 

This function is analytic on the entire complex plane except for a cut on the positive real 
axis in the interval [a 2 , b 2 }. Then one has 

H s (\±ie) =R S (\) Tinps(X) , a 2 < A < b 2 , (28) 



where R S (X) is, from eq.( 



71 

RsW = W ^- 1 (X-- 2 r- 1 - (29) 



Using the standard method of solving the integral equations [IB|, one can show that the 
expression 

Hs{x) = ^x^i u-b 2 x-^ (30) 
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has the correct analytical behavior in f3 < (3 C region (see IJ for more details). Here the 
contour c is a contour encircling the cut [a 2 ,6 2 ] and excluding A. Deforming c to a contour 
around the point A and the contour (a contour at infinity), one finds 



2m V A - a 2 J Cx V ( ~ b 2 A - f 
Inserting ( p9|) in (fTfl), it can be easily shown that p s (A) is (using fl28f)): 



(32) 



At n = 1, where gSCM reduces to SCM, it can be shown that (|32"D is equal 



which is the same as one calculated in |TD| . To find the parameters a and b in eq. (1321) , we note 



that at |A| -> c», eqs.Q27|) and (0) imply #(A) -> 1/A or ^(A - a 2 )/(A - fe 2 )#(A) -> 1/A. 
Therefore we can expand w (A — a 2 )/(A — b 2 )(R s (X) —mp s {\)) and take the coefficients of A 



and 1/A equal to and 1, respectively. It can be see that the coefficient of A is identically 
zero, and the second condition yields to (for (3 < (3 C ) 

^rr E (-1)* (" 7 l ) B p C n ^ p a^b^-^ = 1. (34) 

In n — 1 case, eq.(|3~4"l) reduces to b 2 — a 2 = 4/3, which is one obtained in flQ| . The same 
calculation for the density function p w (X) yields to 

{0 n = 1 

^ TV /(6 2 -A)(A-a 2 )E P , s=0 C p C n ^ 2 X s a^b 2 ^~ s ~ 2 ^ n > 1 (35) 

and the following equations which specify the parameters a and b in j3 > /3 C , 



2(3 2n 



— E C p C n ^ ia 2 % 2 ^-^ = 0, (36) 



p=0 

n 



2(3 2n ~ 



jY.CpCn-pa^b 2 ^ = 1. (37) 

p=0 



Now let us focus on critical point (3 = (3 C in which a = a c and b = b c must satisfy eqs. fl34|) , 
(|36| ) and ( |37D simultanously. It is not difficult to see at /3 = /3 C , except the first term of eqs. 



(134]) and (|37 ) which is equal, all the other terms are not the same. The only unique solution 
of this inconsistency is 

a c = 0. (38) 
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Inserting (|38|) in eqs.(|36D and (|37|) (or (|34l)), results 



n 



V2n-2 

-0, (39) 



2n-l 

c 



2A 



C n 6f = 1. (40) 



2n-l w «"c 
c 



The solution of these equations are 



b c — ► oo, 

Pc — > oo, (41) 

such that eq. (|4"0D holds. Now in the weak-coupling region, where p w ~ 1//5 2 ™ -1 , (3 is always 
greater than (3 C = oo, so p w — > 0. Therefore there is only one regime at d = 0, i.e. the 
strong-coupling regime. 

To check eq. (|22|) , it is easier first to calculate the contribution of Z 0)n to internal energy 
Uq^ and then find Z 0>n from it. If we denote this contribution by U^l, it is equal to (see 
eq.'dl) 



II 



2p 2n 

2n - 1 r b2 



2J , n / a (A-a 2 )>W(A)dA, (42) 
in which pW(A) is given by eq.(B^) and the relation between a and b can be obtained from 



(|34j) . To see the consistency of the results, it is sufficient to check some specific examples. 
In n = 2 , (|^) reduces to 

£/W = _27 ( 6 2_aY, (43) 
0,2 256/? ?V ; V ; 



and ( 34]) reduces to 



therefore 



|(!> 2 - a 2 ) 2 = P\ (44) 



This equation gives Z Q)2 ~ exp( ln/3) , which is the same behavior as eq. (|22|) . It can be 



also shown that for n = 3,4, and 5, these two methods have the same result. As another 
reason for the fact that only the strong-coupling regime exists in d — 0, it can be seen that 
the Z 0jn calculated from p w (X) (eq. (|35|) ) does not coincide with eq . (|22|) . 
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4 Phase structure of gSCM at d = 2 



As mentioned earlier, the gSCM can be solved analytically only in d = and d = 2 cases, 
which the latter is an important case because of its equivalence to gYM 2 (and in special 
YM 2 ). So in this section we want to study the gSCM at d = 2 for V = Tr(AA^) n . 

It is clear from eqs. (|l3|) and flllD that the density function are known in the weak and 
strong regimes at d = 2, as the second term in the right-hand sides of these equations 
vanishes. It can also be shown that at d = 2, p w (0) must be equal to zero at critical point 
(3 = j3 c . In this way, one can obtain the precise value of a, b, and (3 at critical point as 
following M: 



a c = 0, 

2n 

Or 



0c 



2n - 1 
2n 



2n 



n(4n - 3)!! 



(46) 



2 2n (2n - 1) 

To study the phase structure of this model, it is necessary to calculate the internal energies in 
both weak and strong regimes. We know the functional form of the density function in these 
two regimes, but there are two unknown parameters a and b in each of these densities which 
must be obtained from the eqs.([H]) and (fL6|) for weak regime and from eqs.([L7]) and (|i~8|) for 
strong regime. Unfortunately these equations are too complicated to be solved exactly, but 
the crucial point is that as we want to study the phase structure of the model, it is sufficient 
to look at the solutions near the critical point. Therefore we expand the equations (|15|) and 
fll6 ) around a c = and b c = 2n/(2n — 1), and find a w and b w in terms of a = (3 — f3 c up to 



second order. After a lengthy calculation one finds 



An - 3 4n - 3 2 
n = rv — rv -I- 

2(2n-l)/3 c 2n(4n-5)/? c 2 
2n An -I 1 „ 

K = 2n-^l + 2(2n-l)P c a -8^l a + '~ (4?) 

Now inserting p w (z) from eq. (|i~3|) , for a m = 5 n ^ m and d = 2, into eq.(^), we arrive at 
US = - ^IZ-i E B r C p C 2m _ 2 a^- P - q -i I(2n + g _ r )_l + (J__l)i, (48) 

Z7T P P,q,r=0 AU P 

in which 

= ^r(3/2 + t) _ oWs 1 3 5 a 

K J 2r(3 + A;) | + Jfe L 2 ' 2 '2 ' b h K J 

where 2 i*i [a, c, x] is the hypergeometric function. Keeping the terms up to second order of 
a^, eq. fl4"8"D becomes: 

= K{T x bt - l -T 2 b^~ l a w + T 3 bt?- 2 a 2 w + . . .) - 1 + (i- - l)i (50) 
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where 



K 



n{2n-\) 



4n— 1 



^3 = E 

9=0 



Ti =^C 2 „_,_ 2 M( M ,0), 

T 2 = C 2n _ 2 M(n,0,l), 

3 11 

-C 2n _ g _ 4 M(n, g, 0) - -C 2n _ g _ 3 M(n, g, 1) - -C 2n _ g _ 2 M(n, g, 2) 



in which 



M(n, g, r) 



T(3/2 + 2n + g - r) 



T(3 + 2n + q- r) 
In the strong regime, the same calculation leads to 



(51) 
(52) 



An -3 (n-l)(4n-3) , 

-a - _ ,, n J a 2 + 



2{2n - l)p c 2n(An - 5)/3 6 2 



2n 



+ 



An - 1 



-a + 



n — 1 



a z + 



2n - 1 2(2n - 1)/3 C 2n/3, 
for the parameters a and 6, and 

&g = K{T x b? + l -T 2 bf^a s + S^r 2 a 2 s + ...)- 1 + (i- - 1)1 

for internal energy. S3 in eq.(|54]) is 



(53) 



(54) 



^3 = E 

9=0 L 



(~ -n)C 2n _ g _ 4 M(n, g, 0) - ic 2n „ g _ 3 M(n, q, 1) + (f -n)C 2n „ 9 _ 2 M(n, q, 2) 

4 o 



• (55) 



Now subtracting the internal energies in two regions (eqs.(|54|) and (|50D) and using eqs.(47) 
and (|53|) , one obtains 

(56) 



Uii-u!g = f(n)(f3-f3 c r + 



where 



in which 



(2n-l)(4n-3)6f~ 1 
J{U> 32 V ^F/3 C 4 -+ 1 



4n-3 ^ 4n-3 1 

4nTi ; -To Ta >, 

2(4n-5) 2 J' 



(57) 



T 4 = C 2n _ 2 M(n, 0, 2) + C 2n _ 3 M(n, 0, 1) + 2 £ C 2n _ ? _ 4 M(n, g, 0). (58) 

g=0 

Eq.(^) shows that all the V = Tr(AA^)™ gSCM has a third order phase transition near 
the critical point (3 = j3 c . Note that f(n) > for all values of n, as we expected from 
the definitions of strong and weak regimes, and also the value of f(n) for n = 2, 3, and 4 
coincides with those in ref . fllH . 



9 



As a final comment, it is interesting to compare the behavior of the solutions with respect 
to n. The only physical comparable quantity, in the case of third order phase transition, is 
the difference 



In Fig. (1), where we have plotted /3 2 ^(AC) n vs n, one can see a saturating value 2 for this 
difference. 
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